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SUMMABY 



General thin—airf oil theory for a compressi'ble fluid 
is formulated as "bdundary problem, for the velocity poten— 
tialy vrithout recourse to the theory of vortex motion. 

On the "basis of this formulation the integral eq.ua-. 
tion of lifting-surface theory for an incompressible fluid 
is derived with the ehordwise component of the fluid ve— ■ 
locity at the airfoil as the function to be determined^ It 
is shovm how by integration by parts this integral equation 
can be transformed Into the Blot— Savart theorem, A clari- 
fication is gained regarding the use of principal value 
definitions for the integrals which occur. 



The integral equation of lifting— surf ace theory is 
Used as the starting point for the establishment of a theory 
for the nonstationary airfoil which is a generalization of 
lifting— line theory for the stationary airfoil and which 
might be called, "lifting— strip" theory. Explicit expres- 
sions are given for section lift and section moment in terms 
of the circulation function, which for any given wing de- 
flection is to be determined from an integral equation which 
is of the type of the equation of lifting-line theory, 5?he 
results obtained are for airfoils of uniform chord. They 
can be extended to tapered airfoils. One of the main uses 
of the results should.be that they furnish a practical aeans 
for the analysis of the aerodynamic span effect in the problem 
of wing flutter. The range of applicability of "lifting- 
strip" theory is the sasie as that of lifting-line theory so 
that its results may be applied to airfoils with aspect ratios 
as I'oft as three. 
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IKSHODUOTIOK 



ffhe theory of thin airfoils may be characterized "by 
the statement that it is the linear theory, obtained through 
simplifications from an exact formulation of airfoil theory 
for a nonviscous fluid. It is based on the assumptions of 
infinitely thin airfoil sections, infinitesimal angle of 
attack, and infinitesimal camber. In spite of these restric— 
tivo assumptions it is plausible and generally accepted that 
the theory reproduces the characteristic behavior of actual 
wings v/ith finite angle of attack and camber rather accurately 
in many prp.ctical problems, 

first contributions to the theory of thin airfoils are 
Prandtl's theory of the lifting line (reference 1) for the 
determination of the aerodynamic span effect for stationary 
airfoils of not too small aspect ratio, and Munk's two- 
dimensional theory of stationary airfoils (reference 2), 

Shortly after Munk's work investigations were made on 
the two-dimensional linear theory of nonstat ionary motion 
independently by Birnbaum (reference 3) and by Wagner (refer- 
ence 4), Birnbp-um reduces the problem to an integsal equa- 
tion by means of theorems on vortex motion while Wp.gaei- ob- 
tains an integral equation by way of formulating the boundary 
problem for the velocity potential, Birnbaum »s formulation 
(in common with the formulation of more general problems "by 
means of the "acceleration potential" presently to be dis- 
cussed) possesses the disadvantage of leading to an integral 
equation with a considerably more complicated kernel than in 
the Telocity potential formulation of the same problem. Solu- 
tions of this integral equation were obtained by Birnbaum by 
means of numerical methods restricted, ho\*ever, to a far too 
small range of values of the iraportant "reduced frequency" 
parameter which was introduced by him. In Wagner's work 
atte-ition is focuesed on transient problems for which the 
integral equation of the problem is solved^ by series devel- 
opments. Both authors deal only v/ith the motion of .•> rigid 
straight— line profile, • _ 
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Subsequently G-lauert (reference 5), on the "basis of 
Wagner's work, obtains an explicit solution for the cp.sc 
of a simple harmonic motion of a rigid profile. This solu- 
tions depends on certain definite integrals v/hich are func- 
tions of the reduced frequency parameter. The definite 
Integrals are evaluated numerically in Olauert's paper for 
still too small a range of the reduced frequency parameter. 

In 1935 P,nd 1936 Ihoodorsen (reference 6), Cicala 
I'efcrencG 7), Ellonberger, nnd Von Borboly (rofersncos 
8 and 9), and Kusstier (reference 10), independently ?<,nd 
publishing in this chronological order, gave the solur- 
tion of the two-dimensional problem for arbitrary motion 
and deformation of the airfoil and obtained explicit ex- 
pressions for air forces and moments, They also found 
that the definite integrals occurring in Slauert's st)ecial 
and their own general solution could be expressed in terms 
of certain tabulated Bessel functions of the reduced fre- 
quency- par.ameter. 

Modified derivr.tions of these results and applications 
have subsequently been published by Garrick (reference 11), 
Ton Earman and Sears (reference 12), Diotze (reference 13), 
Schwarz (reference 14), S'ihngon (reference 15), and others. 

Approximate solutions for the two-dimensional motion 
of a rigid profile in a slightly compressible fluid have 
been given in 1938 by Fossio (reference 16), 

An account of the work regarding the effect of spc-nv/ise 
variation of the flow (three-dimensional theory) may bo sub- 
divided into two parts, 

G?he first part includes investigations having the pur- 
pose of improving lifting—line theory for stationary motion 
on the basis of various formulations of lifting— surf ace 
theory. Workers in this field have been Blenk (reference 
17), Burgers (reference IS), Von Karman (reference 19), 
Schlichting (reference 20) Tuchs (reference 21.), Bollay 
(reference 22), Wieghardt (reference 23), Kinner (reference 
24), and Erienes (reference 25). 

(The second part includes investigations having the pur- 
pose of obtaining generalizations of lifting— line theory. to 
problems of nonstatlonary motion. Attempts of such general- 
izations on the basis of vortex^-f ilament considerations were 
made independently by Cicala (references 26 and 27) and Von 
Borbely (reference 28), A study of the special case of an 
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infinite airfoil with uniform chord, undergoing "bending 
deformations varying sinusoidally in the direction of the 
span is due to Sears (reference 29), Results "based on the 
thoory of the acceleration potential were given by Zussnor 
(roforence 30), The case of a rigid airfoil with elliptical 
plan form was investigated "by Jones (references 31 and 32), 
IThe author "believes that all these attempts contain assump- 
tions, basically or analytically, which may lead to consid— 
era"ble errors in the determination of the aerodynamic span 
effect, Shis is discussed in the "body of the present paper. 

Lastly, mention is made of two general formulations of 
thin— airfoil theory "by Frandtl (reference 33) and Eussner 
(reference 30) where'by an integral equation is o"btained for 
the distribution of pressure over the airfoil. Since, ac- 
cording to the ecLuatlons of motion for a nonviscous fluid, 
the pressure may be thought of as the potential for the 
acceleration field this approach has become known as the 
acceleration— potential method. An Important point of this 
method is the avoidance of the explicit introduction of the 
trailing surface of velocity discontinuity ("trailing vortex 
sheet"), Zt is felt, however, that nonetheless this method 
possesses serious disadvantages when applied to problems of 
nonstat ionary motion, compared with a method making use of 
the velocity potential. In this connection it is the author's 
opinion that the failure of Birnbaum (reference 3) to obtain 
a complete solution of the two— dimensional problem is largely 
due to the fact that a method was used by him which is iden- 
tical v/ith the two-dimensional form of the acceleration po- 
tential method. 

The present paper is composed of three parts. In part 
I the known fundamental differential eq^uations and bovindary 
conditions of thin airfoil theory are rederived. On the 
basis of these equations the boundary problem for the veloc- 
ity potential for a slightly compressible fluid is formu- 
lated. In this formulation any reference to the theory of 
vortex motion is avoided by means of a simple symmetry con- 
sideration. 

In part II an integral equation is derived for the chord- 
wise conponent of the velocity of the fluid at the airfoil, 
restricting attention to the case of an incompressible fluid. 
This integral equation has the important property that its 
k'ernol is the same for problems of stationary and nonstation— 
ary notion. In contrast to what occurs in the integral equa- 
tion for the acceleration potential. It is shown how the 
integral equation obtained can be transformed, by integration 
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by parts, into a form equivalent to the Blot— Savart theorem. 
It may be mentioned that a clarification is gained regarding 
the use of principal value definitions for the integrals 
occurring in the different forms of the integral equation 
of the problem, 

la part III "the integral equation of lifting-surface 
theory is used to establish a theory for the nonstat ion.-,ry 
airfoil v;hich is a generalization of liftiag-line theory 
for the stationary airfoil and for which the name of lift- 
ing-strip theory is proposed. This nev; theory, which in- 
cludes the known two— dimensional theory for the nonstation- 
ary airfoil as well as three-dimensional lifting-line theory 
for the stationary airfoil as special cases is believed to 
be the first correct theory of this kind. It permits deter- 
mination of the aerodynamic span effect for airfoils of not 
too small aspect ratio in a manner which is a combinntion 
of the known procedures in the two— dimens i onal theory for 
the nonstat i onary airfoil and in the three— dimensional theory 
for the stationary airfoil. Explicit expressions are given 
for section lift and section moment in terms of the circu- 
lation fiinotion \*hich for any given wing deflection is to 
be determined from an integral equation which is of the type 
of the equp.tion of lifting-line theory, Ihe calculations 
may be extended so as to obtain an explicit expression for 
the aileron hinge moment. 

This paper forms part of a project of the Departwont 
of Mathematics at Massachusetts Institute of Technology 
made possible by financial assistance from the National 
Advisory Committee for Aeronautics. 

The author wishes to express his appreciation to Prof, 
H. B. Phillips, head of the Department of Mathematics nt 
M.I.T.jfor relief from part Of his regular duties in connec- 
tion with this work; to Professors H, Reissner, K. Rauscher, 
and f. B. Hildebrand for helpful discussions on the subject 
of the paper; to Prof. Rauscher also for the original sug- 
gestion to study the problem of the aerodynamic span effect, 

LIST OP SYMBOLS 

U velocity vector, giving the difference of fluid 

velocity and velocity of flight 

i?o) velocity of flight ("undisturbed" velocity) 
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U, V 

9 9 


velocity components 




P+ Poo 


pressure 






undisturbed pressure 






density 




Poo 


undisturbed density 




X,Y,Z 


Cartesian coordinates 




T 


time 




P 


implicit eq.uation of the airfoil 


surface : 


H 


explicit equation of the airfoil 
Z = H(X,Y,a!) 


surface : 


a 


velocity of sound in undisturbed 


flow 


Ha 


region occupied by projection of 
on X, Y — plane 


airfoil 



X^ X-coordinate of trailing edge of airfoil region 

x,y,z dimensionless Cartesian coordinates, x a X/b, 
Y = y/b, z = Z/b 

b reference length, for airfoils with uniform chord 

identified with the semi— chord 

t dimensionless time, t s ujlT 

lu reference frequency 

h dimensionless explicit equation of airfoil surface 

h = H/b 

u dimensionless velocity vector, u = V/\.TJ^\ 

p dimensionless pressure, p = P/--p ^ 

k "reduced frequency" parameter, k = bu)/Uoo 

^ Mach's number of the undisturbed flow, ^-Uoe/a 
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x-coordinate of trailing edge of airfoil region 

Idy dimensionless air force associated with strip dY ' 

mdy dimensionless moment of air force associated with 
strip dY 

s 

LdY air force associated with strip dY, L = ^2p^U^ 'b)l 

MdY moment of air force associated with strip dY, 
K = (2p U m 

00 09 

4> velocity potential 

hk amplitude function defined by h = e^* 

«!> defined "by <i> = e^* 4^ 

Pjj. defined by p = e^* pj^ 

x-coordinate of leading edge of airfoil region 

v/ake region, being the semi— infinite strip in the 
x,y— plane extending from the trailing edge in 
the direction of the main flow 

Ey region of x,y— plane which is not part of airfoil 

and wake region 



'k 

it 



circulation function defined by 

""^M^ZjlI^ dx 
ox 

'XI 



2 jjr. 



i t 

Uge chordwise fluid velocity component at airfoil, 



^0 = 



^«t>i]j![x,y, + o) 



bx 

r. defined by T = e^* T. 



w^e normal fluid velocity component at airfoil, 



Zs 0 
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Iji a symbol designating various integrals in the 

course of evaluation or transformation 

defined ty \ = e^* \^ 

mjj. defined "by m » e^* mjc 

s ratio of semi— span and semi—chord for airfoil 

of rectangular plan form 

y* dimensionlesB coordinate defined by y* ss y/s = Y/sb 

~ a sign designating quantities of the two— dimensional 

( section— force) theory 

S a function defined "by equation ( 133 ) 

^ a function defined by equation ( 135 ) 

Eq,^]^ modified Bessel functions of the second kind 

"'^o»'^ii"^2 Bessel functions of the first kind 
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I,- THE BOUNDARY PROBLEM FOR THE VELOCITY POTENTIAL 
IN NONSTATIONARY THIN-AIRPOIL THEORY 



PORMULATXON 07 THE PROBLEM 



iTeglecting finite thickness and finite angle— of-attacfc 
effects, thin— airfoil theory treats iirings as almost flat 
plates, possessing no thickness. Plow of an ideal compres- 
sible fluid is assumed, the uniform velocity of which is 
disturbed by the presence of the airfoil, which is only 
slightly inclined against the direction of the undisturbed 

velocity. The velocity change U caused "by the presence 
of the airfoil is considered small compared with the un^ 

disturbed velocity Ues and the changes in density and 
pressure, p and P are considered small compared with the 
undisturbed density and pressure p„ and Pq,, On the basis 

of these assumptions a linearized form of the problem is 
obtained. 

Before linearization the differential equations of 
the problem are the Euler equation 



li!jL!£^.(i?+$^).grad (U+U«)«-ii:^-i!^^l (l) 
and the equation of continuity 

— ^T^ + C(p+p^) (U+Tf„)3 = 0 (2) 



As boundary condition it is prescribed that on the 
surface of the airfoil the normal component of the fluid 
velocity-- equals the normal component of the velocity of 
the element of the airfoil with which it is in contact. 
If P(X,Y,Z,T) - 0 is the equation of the surface repre- 
senting the airfoil the boundary condition has the form, 

P(X,Y,Z,T) = 0; ~ + (U+T7„).grad P = 0 (s) 

hT 
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Formulation of the conditions at infinity and of the con- 
ditions along the edge of the airfoil surface is postponed 
until after the problem has been linearized,^ 

Ihe undisturbed velocity may be taken as parallel 

to the X-axis, so that 

GO CO 

where "i" represents a unit vector in the X-direction, 

Ihe airfoil surface is assumed to lie very nearly in 
the Xjl-plane (fig, l) and its equation may be written as 

P s Z - H(X,Y,T) = 0 (5) 

The asBumption of small disturbances is equivalent 
to the following order of magnitude relations, for the 
.velocity changes, 

t 

IU| Ueo, U« (6) 

For the slope of the wing surface, 

— <*: 1, — « 1 (7) 
bX ' bY 



The condition of small density change is 



(8) 



permitting the linear pressure change density change re- 
lation 



^a?he stationary problem (b/bT^O) in its general form 
has boon discussed by R, von Mises in a paper given before 
the American Mathematical Society in April 1942, 
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where 



is the square of the velocity of sound of the ■undisturbed 
flow. 

Introducing equations (6) to'(lO) into the differential 
equations (l) and (2) and neglecting terms small of higher 
order leads to the following linear equations, 

-I- + tt.grad U = - grad P (ll) 

+ p div TJ + div (pU«,) = 0 (12) 

vfhile the "botuidary condition (3) becomes with grad Z=k 
where ^ stands for a unit vector in the Z— direction, 

_ 15 + U.k - U«>«grad H = 0 (l3) 

In viev; of the fact that the airfoil surface lies very 
nearly in the X,Y — plane it is permissible to satisfy 
this condition instead of at the surface itself at the 
projection of the surface onto the X,Y— plane. Denoting by 
the region occupied by the projection of the airfoil 

surface, which will henceforth be called the airfoil region, 

and introducing for grad H its value i ~f + j" ~ and 

oX 

writing V/ for the Z-component of the velocity, the 
boundary condition takes on its final form 

Z = 0; X,Y in R^: W = II + Uoo 11 ( 14) 

The region E^^ is to be considered as the limit of 
closed surfaces surrounding Sa, the boundary c ondition 
(l4) holding on both sides, Z = ± 0, 



Concerning the shape of the region R„ the restriction 
is made that straight liijes having the direction of 
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intersect its boundary either in two points, are tangent to 
it, or are entirely outside. 

In addition to the differential equations (ll) and (12) 
and the "boundary condition (14) the fundamental condition is 
imposed that along the trailing edge of the airfoil tho ve- 
locity remains finite (Eutta-Joukowsky condition) 

X^j (y) a Oj U finite (15) 

It should "be remembered that this condition is motivated 
in the follov.'ing way. Experiments show that for airfoils of 
finite thickness with rounded leading edge and sharp trailing 
edge the effect of viscosity manifests itself in such a way 
that the flow pattern is cLuickly developed into very nearly 
the same as that for an ideal fluid with the condition im- 
posed that the velocity remains finite at the trailing edi^-e, 
which is the only place where it could mathematically become 
infinite. Conditions for the existence of this close connec- 
tion between viscous and ideal fluid flow theory involvo 
limits on the thickness-chord ratio of the -airfoil sections 
and on the magnitude and direction of the undisturbed veloc- 
ity U(» which are satisfied for conventional airfoils with 
angle of attack below the stalling angle and velocity of 
fUsht not too close to the velocity of sound. 

Ho condition is imposed in thin-airfoil theory on the 
velocity at the leading edge, A sharp leading edge is con- 
sidered as the limit of rounded edges for which the velocity 
at the leading edge becomee, in general, infinite. Ho the 
extent that this happens, linearized thin-airfoil theory 
must be considered inconsistent. The excuse for permitting 
such inconsistency is furnished by the fact that it is re- 
stricted to a zone of small width adjacent to the leading 
edge, This makes plausible that the effect is insignificant 
so far as it concerns the calculation of the resultant forces 
and moments which the flowing fluid exerts on the airfoil, 
(In the two-dimensional stationary theory this has been con- 
firmed by comparing the results of the linear theory i/ith 
the known results of exact theory.) 

Conditions at infinity are undisturbed flow far in front 
of the airfoil — that is, 

Xpr-m; U=0, P=0 (16) 
while, as regards conditions far behind the airfoil (Xs+os) 
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the omission of the viscous teriss in the Euler equations is", 
in general, responsxljle for a persistence of the disturbing 
influence of the airfoil, 

BiMEIirSIOHLESS FOBM OF THE S(iirA!ZI OUS OF SHE PROBLEM 



Before proceeding further it le convenient to make the 
system of equations (ll) to (16) dimens i onles s. Dimension- 
less coordinates are introduced "by putting 



X y = I, a » I (17) 



where "b^ is a reference length which in the two— dimensional 
theory vrill he identified v/ith the semi-chord of the airfoil, 

Simensi onlees time is introduced by putting 

t = u) T (18) 

where u) is a reference frequency \*hich in the case of 
harmonic oscillations will be Identified with the frequency 
of oscillation, 

A dimensionless camber surface equation is introduced 
by putting 

h=^ (19) 

Dimensionless velocity and pressure changes are intro- 
duced by putting 

— > 

- - ^ (20) 

p = (21) 

When equations (17) to (21) are introduced into the 
differential equations (ll) and (12) and into the bound- 
ary conditions (l4), it is found that it is convenient to 
define the following two dimensionless parameters, 



:TACA tit No, 946 



14 



k = Ji5 (23) 

P = (23) 

a 

Uhe "reduced frequency" parameter k is of "basic importance 
in the theory of nonstat ionary motion. The parameter p 
represents Mach's mimber of the undisturTied flow, the ratio 
of the velocity of the undistur'bed flow and of the velocity 
of sound in the undistur'bed flow. 

She differential equations and "boundary conditions have 
now the following form, if use is made of the relation, 

k||4.|?=-i gradp (24) 

1 ^« f'k ^ + div ^ = 0 (25) 

2 V at ax y 

z = 0; x,y in Sa; « = ^ |^ + ^26) 
z = 0; X = X|.(y); finite (27) 
X = -00 ; u 5= 0, p = 0 (28) 



Che analytical part of thin— airfoil theory comgists 
in solving this system of equations (24) to (28) so that 
e:cpressious may he obtained for forces and moments exerted 
by the fluid on specified portions of the airfoil, IDhese 
forces and moments manifest themselves through discontinu- 
ities of pressure v/ithin . the airfoil region H^^, A conven- 
ient notation is to write p^ for the pressure on the side 
2 = + 0 of Sg^ and p^ on the side z = — 0 of IL^, 

Denoting by Idy and mdy dimensi onless forms of 

forces and moments associated with strips dy of the a-ir- 
foil, a convenient way of writing forces and moments pei- 
unit of span length is 
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I (y; Xj^, Xg) = i / [p^ (x,y)- p^(x,y)3 dx (29) 
1 r 

3Cp, x^,X3) = ^/ (x-X^) .Cp;(x,y) - p^(x,y)] dx (30) 

where from now on the variables y, , Xg In \ and m 
will not Tse written unless necessary to avoid ambiguities. 

!I!he relations between dimensionless forces and moraents 
and the corresponding actual forces and moments L and H 
are established by means of equations (31) and (17), 



1 = C2p^ Ujb)l ■ (31) 

K = (2p„ U„^b^)m (32) 

EEDirOSION Off THE BOUNDARY PROBLEM TO A 
PROBLEM FOR THE HALI" SPACE 



So far the problem has been formulated as boundary 
probleu for the exterior of an infinitely thin closed sur- 
face. It may be observed, however, that ' in e quat ions ( 24) 
to (28) of the problem no distinction occurs between the 
portion of the fluid "above" the airfoil surface and the 
portion of the fluid "belov/" the airfoil surface, This 
indicates that the flow must possess properties of symmetry 
and antisymmetry with respect to the x,y— plane, in the 
sense that the components of the velocity vector u and 
the pressure p are either even or odd functions of the 
z— coordinate, The boundary equation (26) which is to be 
und or st ood as holding for z s diO indicates that the ve^ 
looity component w is an even function of z. Prom the 
z—componont equation of equation (24) it follows then that 
also bp/bz is an even function of z and consequently p 
an odd function of z, Prom the remaining component equa- 
tions of equation (24) it follows then that also the volocity 
components u and v are odd functions of z. 

The fact that the pressure p is an odd function of z 
in conjunction with the fact that p is continuous within 
the fluid indicates that p vanishes over the part of the 
x,y-rplane outside the region Ra, 
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z = 0; x,y outside Ea s p = 0 (33) 

OoHbining equation (33) with equations (26) to (28) there 
is obtained a system of conditions for the entire x^y— plane 
which may then he considered as the boundary of one of the 
half spaces, say the half space z > 0, 

lETTRODUGriOir OF VEIOCITY POTENTIAL 



Further treatment of the problem is carried out here 

in terms of a velocity potential <j>, the existence of v/hich 

follows from the fact that equation (24) implies 

k JL + A )rot u*= 0 (34) 
^ St Sx / 



Since rot 0 in the region of undisturbed flow, it 

f'ollovrs from equation (34) that throughout, the interior 
of the half spaces z > 0 and z < 0, 

rot 1? = 0 (35) 

Prom this it is concluded 

Ti = grad * (36) 

Introducing equation (36) into the continuity equation (25) 
leads to 

2 \ at bx/ 

Introducing equation (36) into the equation of motion (24) 
leads to 

grad (^k ^ + - - grad p (38) 

which is equivalent to 
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1 , b4> , , 

with an additive arliitrary function of t incorporated in 

Introducing eq.uation (39) into eq_uation (37) furnishes 
the differential equation for 4>, 

a 

p= (k± + ± ) * = 0 (40) 
\ ot ox ^ 

In terms of <b the boundary conditions (26) to (28) 
and (33) assume the follovring form: 



z = 05 x,y inside = k + — (41) 

bt ax 



z = 0; X = x^(y):: grad 4) finite (42) 

z=0; x=-«';'(|)-0 (43) 

z = 0; x,y outside Ha_ ? k~S~+--- = 0 (44) 

ot bx 



^0 simplify the work from here on the case of simple 
harmonic motion is assumed hy putting 

li(x,y,t) = hj^ (x.y) e^* (45) 

Because of the linearity of the problem more general solu- 
tions may he obtained from the solution for this case by 
means of superposition, for instance, solutions for tran- 
sient problems through the means of Four ier— Laplace inte- 
grals, She functions h^ may be complex themselves so as 
to account for phase differences between different hjj., 

With h given by equation (45a) the potential <j> 
will be of the form 
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4> = 4>jj. z)e^* (46) 

whore again ^y. may "be a complex function. 

Substituting ^ of eq.uation (46) in eq.uation (40) 

*k - (i k + A^c>j, = 0 (47) 

SuTjstituting (}> of eq^uation (46) in equation (39), 

1 1 it / ^ 5<bi>. it , 

The "boundary conditions (4l) to (44) become 

z = 0; x,y inside Ba ' : = ikhfe + (49) 

dx 

z = 0; X = x^(y)-: grad finite (50) 

z = 0; X = - CQ • : 4»j^ = 0 (51) 

z = 0; x,y outside R ■; il«l>]j. + = 0 (52) 

bx 



It is imporisant to no.te that equation (52) may be 
integrated to 

(=c,y,0) = c (y) e-^^^ (53) 

when X and y are out side "the region Rg^, Using equa- 
tion (51) it is seen that c vanishes along lines y = 
constaat vjhich do not pass through the regi'on and that 
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c vanishes along the part of lines y = constant passing 
through reaching from x = — «» to the leading edge 

coordinate , Hence if in addition to the airfoil region 
Hg^ there is defined (l) a wake region B.^ Tseing the seni- 
infinite strip extending from the trailing edge in the di- 
rection of the main flow and (2) a region denoting the 

remaining part of the x,y— plane, the boundary condition (52) 
may "be formulated in the following form (see also fig. 2) 



2 = 0; x,y 



z = 0; x,y 

where Ty,- .stands for the value of ^15. at the trailing 
edge 



and v;.here, since 4>i5.(x,y, + 0) = - <t> ( x, y, - 0) , T = Tjj-e^"^ 
is the value of- what is known as the circulation. 

The problem is now to determine the solution of the 
differential equation (47) for the half space z > 0 
suhjcct to the "boundary conditions (49), (50), (54), and 
(55), The notable feature of this mixed boundary problou 
is that it contains the undetermined function Tjj. which 

has to be found in such a form that the finiteness condition 
(50) is satisfied. 



in . Rr ' *k 0 (^4) 

« . A Ir / N -ikC=E-xt(y)] . . 
in R„:J •• *ic » "5 rk(y)e (55) 
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II,- IHB ISTSaaAL a(4UATI0N or LiyTING-SUSffACE 
THEORY FOE AN INO OMPRBSS I3LE PIUID 



Equation (47) indicates that for an incompressible 
fluid c{5jj and also 'b^/'b'x. are harmonic functions. The 
function^ h^/'dx. may be represented in the interior of 
the half space z > 0 by means of its values at the bound- 
ary, z s 0 which may be denoted by Ug, Dispensing in v;hat 

follov/s with the subscript k this representation is 



il* »= i rr ^qCI. n) z dt^ ^ 



37i 



(57) 



The kernel 



K = _L ! g-T^g (58) 

2" TCx- |)% (y-n)^+ z^j 

may "be obtained either by way of introducing spherical 
coordinates about the point (|»Ti,0) in the differential 
equation for h^/hx. whence a solution possessing the proper 
singularity is found by separation of variables, or by means 
of a Jourier integral solution for b4>/^x for the case that 
Uq vanishes everywhere except over an infinitesimal .area 
dldti v/here it has the value l/d|dTi, 

To obtain on the left hand side of equation (57) the 
correct boundary values Uq when z = 0 the integral aust 

be defined properly. One such definition is obtained in the 
following way. Write equation (57) in the form 



^ -If ^ If 



where Ej^ is a small rectangular region surrounding the 
point I = X, Ti as y, and Ej^ is the remaining region of 
integration^ It is found that the correct boundary value 
Uq is obtained fyom eq.uation (59) when the sides of 
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and the value of z are Simultaneously made small in such 
a wa" that, the value of z tends more rapidly to zero than 
the smaller of the sides of the rectangle H^, This result 
holds independently of the location of the point x,y with- 
in Hj^ so that no reason to prefer the principal value of 
the integral — which corresponds to symmetrical location of 
x,y within . S.^ — appears at this stage. The motivation 
for working with principal values appears at . a later stage, 

The integral equation for u is obtained by means of 
the condition that over part of tfie boundary the values of 
the normal velocity b(>/c>z are prescribed, The value of 
h^/hz is determined by first taking 



dxdz 



( 



r--37rK^<in (6 0) 



where it is legitimate to differentiate inside the sign of 
integration as long as z > 0. The value of b<J>/bz is 
obtained by integrating equation (60) with respect to . 

This integration is facilitated if use is made of the follow- 
ing identity, 

.35 



02 



Substituting equation (61) in equation (60) and integrating 
from X = —00, where &4>/ &z = 0, there follows 



hz 2tt J J ° 1 0 . P^®^/' 
I (xt_ I) + (y-T^ 



) +z y /_oD 



by 



I 



(y-Ti) dxi 



(xi- i)^+(y-Ti)* + 



rr37a 



7 



1 



dldn 



and after evaluation of the inner integral, 
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r 



hz 2-nJ J °)f. , « 

^ ... i (x-l) +(y-Ti) +z Y 



J 



i S 2 

( (y-n) +z 



S 2 



1 



IcLtcLn (63) 



In this integral the region of integration consists of the 
airfoil region Eg, and the wake region R. 

tion (54) indicates that u. vanishes in B; 



^ only, as eq.uar-. 



To obtain the integral equation for u^^ the boundary 

conditions (49), (54), and (55) have to "be used and the 
limit of the integral (62) as z tends to zero has to be 
taken. In principle the simplest procedure would be to 
evaluate the integral first for non-'^vanishing z and then 
make z = 0 in the integrated expression. However, inas- 
much as the integrand in equation (62) becomes of simpler 
form *\rhen in. it z is made equal to zero it is desirable 
to define the process of integration in equation (62), or 
in an equation derived from it, in such a way that the 
correct result is obtained if first z is made equal to 
zero and then the integration is carried out. This, as will 
be shown, is possible without transformation as regards the 
first term in equation (62) while the second term has to be 
brought into a different form, in order to avoid singular** 
ities of too high order. 



In this transformation two different cases are dis- 
tinguished; (1) the case for which the leading edge is 
straight and perpendicular to the direction of the main- 
flow, (2) the case where this is not so. In the first 
case the second term may be integrated by parts with re- 
spect to r\ and since in view of equation (54) u^^ 
vanishes at both ends of the ti— integrat ion interval 
there follows from equation (62) 
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1 

2tt 




/ / S'^o t: 



( 



>d|dTi (63) 



In the second case the" transformation of equation (62) 
involves first integration Tjy parts with respect to |, 



CO "» 



d| 



Prom equation (54) follows that in the integrated part -<j>o 

vanishes v;hen t = ^\ v;hile for | the factor of 

vanishes. Thus instead of equation (62) there may "be 
writt en. 



r 

I ".(x-l) + (y- Ti) +z j- 



oy 



y - Ti 



X — 



(y-.0«+z« .(x-.|)%(y_^)%,V 



i/a 



+ 1 



In this expression the second term may "be integrated "b.^' 
parts uith respect to ti and, writing b4>Q/'br\ - v^ 

and observing that the integrated part vanishes at the 
limits, there follows 
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b4> 



1 

Sir 




r 



3/a 



IdldT, (64) 



It may be noted that an equation equivalent to equation 

(64) could also have "been oT:tained by applying the Biot— 

Savart theorem with a vortex sheet occupying the regions 

E and R. in the x.y— plane, 
aw 

Hhe integral equation of the problem is now obtained 
by substituting either in equation (63) or in equation ,(64) 
the boundary conditions (49) and (55), There follows, from 
equation (63), 



x,y in : ^ 



B 



bu. 



L 



a 



(x- I) 



u. 



"sTa 



■^^t) +(y-T\) +z J- 
X - I 



5^ (y-n) +z 



/(=^|)%(y-n)^=s^r^''^ 



+ 1 



d|dTl 



r(Ti) 



-ik[|-Xt(T^)] 

e (x-l) 



2 S 



-ik| d 



v. 

l^r(n)e *j (y-Ti) / ^ _ I 



r 3 /a 



(y-n)^+ 2^ 



{(s>.|)2+(y-.Ti)^+z^-* 



+'1 



d|dn 



.(65) 
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Equation. ($5) is valid when the leading edge of the airfoil 
is a straight line parallel to the axis. ^rom equation 
(64) follows, without restriction concerning the leading 
edge curve, 



X,:/ in R i 



r 



jj^a 



u 



+ V. 



y- n 



d|dn 



r 



I// 



- ikr(Ti) 



,ik(|-Xt) 

e ( X - 1) 



' 3 /a 

2. 2 L 



+ Afre^^^t 



e (y- Ti) I 



^^.(x-|)%(y-t^)%z^} 



3 /S 



(66) 



J 



In "both equations T is defined, according to equation (56), 
ty 



|r(n) =f % (I , Ti) d| 



(67) 



According to equation (50), the additional condition is 
imposed that 



•^o n) finite 



(63) 



Equation (6.5) has the advantage that in it only Uq occurs, 
"but the disadvantage that its applicability is restricted 
to airfoils with leading edge straight and perpendicular 
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to the direction of the main flow, EcLuation (66) is free 
of this restriction, "but contains the two unknowns Uq and 
Vq (vfhich, however, are Tioth derivatives of the same func- 
tion ^o)« developments of part 111 of this paper are 
"based on equation (65) but it is certain that equivalent 
results can be obtained on the basis of equation (66). 

So far as the integral equations (65) and (66) are 
concerned it is possible to put z directly equal to zero 
in some of the wake integrals as the variable x is alv;ays 
exterior to 2^^, Ihis is not the case with the integrals 
over the airfoil region where the integrand becomes infinite 
when z = 0, | = x, y = n* order of this infinity is 

most easily recognized if cylindrical coordinates about 
the point x,y are introduced, that is, 

I — X =! p cos 6, n — y = p sin 8, d|dTi= pdpd9 

whence, for instance, 

(x_ I) d^dn - cos 6 p'^dpde 

((:^|)«+(y-.0'^+z«| [P ^ ) 



On the basis of this formula the integrals over and 
Rg^ could be considered as composed of the following t\io 
regions? (l) a small circle with center p = 0, (2) the 
region minus. the small circular region. It is possible to 

shovr that lim / / has the value which would also be ob— 

J J 

tai.ied by excluding from the region of integration the small 
circle, assuming in the integral over the remaining region 
z equal to zero and by finally letting the radius of t'he 
excluded circle go to zero, 

Por the applications it is, however, more convenient 
to subdivide the region B.^ in a different way. This nay 

be explained in detail for one of the integrals occurring 
in equations (65) and (66): namely, 
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The integral Ig may "be written, as 

h-fj* If 

where H.^ is a small rectangular region surrounding the point 
I = x,7 = T\ in the |,n-plane and Eg the remainder of 
the airfoil region Eg.. la equation (70) Rj^ may he nadc 
so small that in it the function u^ changes very little, 
so that 

r r n n (^c-l) d|dTi 

/ / = // (71) 

\fith Ej^ thus determined make z so small that 



u^ce.n ) (x-n dgdT^ _ Up(i,n)(3^naidTi ^ 
Eg |(^e)^+(y-Ti)=+=!«} } 



3?hus 

n ^ (X~|) d|dTl 
= CTio(35:,y) + 53^3 / / 3/k 

^1 "l^^^-^^ +(y-Tj) +z f 



"^•^ ^(x-e)%(y-T,)H 
Ep V. J 
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Since in this definition z approaches zero automat icall;"" 
as the size of the rectangle Rj^ shrinks, there may "be 

wr itten 



:o = lim l^ = iini 
z -♦0 Ej^-*o 



= Uo(x,y) 1 i m / / 

S.->o J o i 



r r ^o<l»Ti) (x-|)d|dti 



+ 1 i 
E 

S 



The most convenient form of Iq is obtained if ther 

shape of the rectangle R^^ is chosen such that the first 
intejrral in equation (74) vanishes. This is fhe case when 
the rectangle B.^ is symmetrical about the line '| = x, 

since the integrand is an odd function of | -x. This is 
the scae as saying that the gap in the second integral is 
symmetrical so far as the |— integration is concerned, 
which is equivalent to saying that the principal value 
of the integral i's taken in carrying out the integration 
with respect to | , This nay he indicated "by writing 



lo 



r / ^o( i. Ti) (x ^1) 

R i (^-^) +(y-^) f 

In the same way it can be shown that the simplest form of 

the remaining singular integrals in equations (65) and (6S) 
is obtained by making the rectangle R^^ symmetrical about 
the line y = n — that is, 

^ r 5uo y - n / ^ - I \ ^ 

- f f ^° -i- ('-^—-iLllS . iLld, (,6a) 



1 i m . 
2— »-o J 
R 
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and 



1 i m 
z— »-o 



R 



a 



3/3 



As in one integral gap symmetry in the |-.direct ion and ^ 
in the other integral gap symmetry in the n— direction is 
desirable, it is most convenient to have in both integrals 
gap symiiietry in hoth directions. With this understanding 
the integral equations (65) and (66) may be written 



x,y in Ra J 



Tio^^.Ti) (x-l) 



w = — 



J 



OUf 



c5n y-Ti 



ik 



|(3«>-l) +(y-T^) j 



+ 1 



didti 



4in 




-ik(|-Xt(Tl)) 

r(Tl)e (s^l) 



(^l)^+(y-Ti)^j 



3/8 



y 



{(x-e)%Cy>,)^}^'^ 



+ 1 



1 



dldTi (77) 
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or 



in H, 



3tt 




< 



J 




-ik(l-xt) 
-ikr(T\)e (x-l) 



e (y- 



-1 



(-n^(y-n)^]''' j 



(78) 



la equations (77) and (VS) P is defined, as before, by- 
equation (£7) and Uq is, as before, subject to the con- 
dition of finiteness of equation (68). 

In what follows equation (77) will be specialized to 
the following cases (leaving deductions from equation (7S) 
for future work): 

( 1) Ihe two— dimensional theory, where no new results 
are obtained; 

(2) iCho stationary theory for airfoils of rectniiguliar 
plan form from which in a manner similn.r to that used by 
Burgorc (reference I's) there is obtained Prandtl's cqur-tion 

of lifting-line theory, Also obtained are expressions for the 
Bpanvrise variations of total moment and aileron hinge moment 
which are equivalent to earlier, apparently not vrell knov;n, 
results of &lauert. 
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(3) The nonstationary theory for airfoils of rec- 
tan^ulc-r plan form. Here a thory is obtained for the 
aerodvnamic span effect which is a generalization of the 
Prandtl theory for the stationary airfoil and is considered 
just as reliable for the nonstationary airfoil as lifting- 
line theory for the stationary airfoil. 

Before proceeding with this program expressions may be 
obtained for lift and moments in terms of the function uo, 
which is the function to be determined from the form of the 
integral eq.uatioa of lifting-surface theory as given in this 
paper. Substituting in equations (29) and (30) the value of 
p froiA eq.uation (48) and taking into account that Pi = - 
there follows first 

&*k \ it 



Xg 



e" dx (79) 



==1 



ni(r;Xo.x,,x,) = y (:^Xo) (ik*^ + ■^)e'* dx (80) 
Ifriting . 

V = Ij^ e" m = mjs: e^* (Sl) 



there follows 



= y (ifc<J^ + Ujj) dx (82) 



3=1 



ana 

X 



(x-Xo)(ikc!>k + Uo) dx (83) 



By integration by parts equations (82) and (83) may be trans- 
formed into expressions depending only on Uq, 

3=1 3=2 «8 

lk(x^,::c3) = ik ^(xg-x^) y uodx+y^ (xg->x)uodxj+y Uodx (34) 
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^ ^X, 



+ 



y^'cx^-x) (fslf-xo) dxj+ y f:^Xo) u^dx (85) 



JroE equations (84) and (85) there follows in par- 
ticular for 

(3.) the section lift (x^^ = xV, x^ = x^) 

Ijj. = ik / (x^-x)Uq dx + / dx (86) 

^xi -^x^ 

(2) the section m oment about Xg (x^^ = =^1 1 = ^t) 

^t 

infcU^) = ik r (xt-x) (^-,xju„dx+ y (x-x^)uodx.(a7) 

(3) the aileron hin ge mom ent (Xj^ = Xq=c, Xg = x^ ) 

c 

dx 



aj, (c,c) = ik |i (xt-.c)»y u^ 



i=t 

+ 



+ 



(x-c) Uq dx (38) 



c 

It me.y be noted that eq.uations (86) and (87) coincide v/ith 
expressions previously used Tjy Olauert (reference 5). 
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III.- LirilNa-STEIP THEORY FOR THE UONSTATIOUARY MOTIOiT 

oi" A3sr AiaroiL of finite span 



In what follows a theory is developed for the airfoil 
of finite span aulajected to nonstationary motion which may 
be considered as a generalisation of lifting-line theory for 
the stationary airfoil. No use is made in this development 
of the vortex filament motion. The starting point of the 
developments is the integral equation of lift ing- surf ace 
theory in the form of equation (77). In this integral equa- 
tion simplifying assumptions are introduced of which it is 
apparent that they are of the same nature for the stationary 
and for the nonstationary airfoil. Thus, the ranie;e of valid - 
ity of the theory -put forward here ooincideB with the range 
of validity of lifting-line theory for the stationary airfoil . 

While derivation of the results for the wing in non- 
uniform motion depends on the same order of magnitude rela- 
tions regarding aspect ratios as the derivation of the re- 
sults for the wing in uniform motion, the steps involved in 
the solution of the uniform-motion case are naturally of a 
much simpler nature than the steps involved in the solution 
of the nonunif orm-motion case. 

The results obtained here consist in explicit expres- 
sions giving lift and moment intensity at every section of 
the span for any deflection of the wing in terms of the cir- 
culation function which has to he determined from an integral 
equation of the nature of the lifting-line equation. If the 
assumption of two-dimensionality is introduced into the re- 
sult s" they reduce exactly to the known results of the two- 
dimensional theory. 

Ho expressions are as yet given for aileron hinge mo- 
ments. Such expressions may, however, "be obtained from the 
present results. 

Also, airfoils of rectangular plan form only have for 

the time being been considered, for the sake of perspicacity. 

It is certain that equivalent results can be obtained for 
tapered airfoils. 
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JOR AH AlRrOlL Off SEOSAH&UlAS PLAN TOHM 



taking as airfoil region R^ the rectangle bounded "by 

the lines x « ±1 and y » ±b, the basic integral equation 
(77) becomes 



_1_ 
2tt 




uoCe.>'n)(x " o 

i (x - I) + (y - Tl) ^ 



9*0 y - Tl 



I 



^ - e 

<;(x . |)« H. (y . T1)H''' 



+ 1 



dldTl 



ik 8 
Ike^^ r r -iki 

4lT 



(y - Tl) 



a* 3 /a 



an 



r 



- e 



y - Tl - |)^ + (y . Tl)^| Jj 



with, according to equations (67) and (68) 

n) d| 
Uo (1,T)) finite 



1 r (Ti) Uo (I. 

-1 



(89) 



(90) 
(91) 



Before considering the three-dimensional problem of the 
nonstati onary airfoil this equation Is specialized for the 
two-dimensional case and for the three-»dimenaional stationary 
case and some results are established pertaining to these 
cases which it is convenient to make use of later on. 



j^ACA 'STS No. 946 



35 



THE IHTEOHAI BCitTATlOir Of IWO-DIHSlHSIONAti THBORY 

The airfoil region has now the form of the Infinite 
strip Ixj < 1 and the aseumption of two-diaensional flow 
is expressed "by 

Vf^ist^y) = Wo(x), Uod, 1) « tt^d), r(Tl) constant (93) 

iSquatlon (89) reduces to 



- - ^ f Trfrr^ '''' 

-x-» iix « |) + (y - tl) > 

iL - £)' + (r-n)'} 

•- 

fChe integration with respect to T) can "be carried out, leav. 
ing 

where T is given in terms of Uq, "by equation (90) and 
where the finiteness condition equation (9l) has to "be ob- 
served. Equations (94), (90), and (9l) can be solved ex- 
plicitly for Uq, the result may be substituted in equa- 
tions (86) to (88) for lift and moments to obtain the well- 
known Theodor sen-Gioala-Ellenberger-K{tssner results. These 
same results will appear as special eases of the new theory 
taking account of the aerodynamic span effect. 

THS XNTS&HAL EQUAO^ZOir OF TEBSS-DZMSUSIOIirAIi S7ATI0HAR7 TH20BY 



Vfith 



k * 0 
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equation. (89) reduces to 



3 /a 



3uo X 
971 y - T) 



{ 



X - 



(x - D + (y - T)) 



21 .i/a 



+ 1 



dIdTl (95) 



while equattone (90) and (91) rcnaln unchangedT 



Approximate solutions of equation. (95) for airfoils 
without eam"ber for which = a = constant have heer. oTd- 

tained "by Blank (reference 17) and Wieghardt (reference 23) 
with the alia of supplementing Prendtl's lifting-line theory. 
The Inverse prolilesi, to ohtaln convenient expressions for 
Vq whon Uq is given, has heen dealt with, using Fourier 
integrals, "by Von Karman (reference 19) and Juchs (reference 
21). 

It appears that the task of obtaining quantitative so- 
lutions of equation (95) giving reliable corrections for the 
results of lifting-line theory is Of considerable difficulty 
and requires work going beyond what has been accomplished by 
Blenk and Wieghardt. 

Lifting-line theory may be obtained from equation (95), 
substantially according to Burgers (reference 18), in the 
following manner. 

Substitute in equation (95) as new variables 



* = 2 

s 



(96) 



v'hich changes equation (95) into. 



, , s rr"{ u,(|,Ti*)(x- I) 

Wo(3c.y) = - ^ of ' ^ a/g 

2" i(x - 1)^ + s^(y* - Ti*)=;>.. 



1 9u. 



s 5T\* y* - Tl' 



/(x ^ + e^(y* - T|*)y^ 



+ 1 



d|dTl (97) 
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For sufficiently large s, practically for 

a > 3 (98) 

the terms 



^ (x - + 6®(y* - Tl*)^!^ 



have a steop maximum for . y* «= T|* so that, roughly, the 
main oontriljution to the value of the integrals containing 
these terms comes from the immediate neigh'borhood of the line 
y* a Tl*. In this noighhorhood the values of vlq and - ^o 

are thought to change sufficiently slowly to permit replac- 
ing their actual values "by their values at y* «* T|* , If this 
approximation is accepted oc^uation (97) becomes 



s-A// r U o(§.y*)(x ) 

L 

1 Suq 1 .X - I 

8* 3y» y* - n* >;(x - V + 8*(y* - nO^V^^® 



1^ 52i2 ^ 
8« an* y* - Tl* 



d|dTl* (99) 



In the first two integrals the integration with rospect 
to T|* may be carried out explicitly. Because the main 
contribution to the value of the integrals oomos from the 
immediate neighborhood of y* = Tj* the error introduced by 
Integrating from -oo to -fm Instead of from -1 to •(•1 is 
neglected, (This latter approximation evidently ceases to 
be good in the immediate vicinity of the tip sections 
y* « ±1, and could not be made were it not for the fact 
that uo turns out to be small near the tip sections.) With 
the following values of the relevant two integrals, 



dT)* - 2 

— ■ -r/s- - 



^09 



|(x - |)^ + s^(y* - Tl*)*j. 8(x - if 

an* " 

(y* - Ti*)^<x - + s^(y* - 
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the Integral et^uatlon (99) le reduced to 

° TTj, X - I ^ 21T8 J J 9T1* y* - Tl* 

-1-1 

and if in the second term use is made of equation (90) there 
follows as final fQrn of the approximate integral equation 
of the stationary rectangular lifting surface. 



1 



» - i f ioiLiIi i{ , J_ / iL (101) 

TT J X " t 4ttbw' dTl* y* - X\* 

It may "be seen that the second term in eciuatlon do/) gives 
the finite-span, correction^ v;hile neglecting the second term 
is equivalent to assuming two-dimensional flow at every sec- 
tion. 

So determine the functions Uo and T use is made of 
a known inversion formula which is to be considered as a re** 
suit of two-dimensional potential theory. (See, for instance, 
Sohngen, reference 34.) She inversion formula states that 
to the relation 

1 

g(x) = - i/-liLl 4| finite (lOS) 

IT J X - 6 

there corresponds the following inverse relation expressing 
f in terms of g, 

fU) /IZEf (103) 

Tt V 1 + X J v'l-l x-6 

-1 

Applying equations (lOS) and (103) to equation (101) 
in order to solve for uoi 



L 

1 r 6.T dTl* 1 A I 
4rreJ ^ dTl* y* - Tl*(x -> i 
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and with 

1 



y-ym' _iL_ . (104) 



5 X - e 



there follows 



1 ~» 



, A./^ ^ (105) 



—1 

Sq.uation (105) "becomes, Integration, an equation for 

1 



ttJ V' 1 + x)J V 1 - I X - I 
-1 J 

It is plausl'ble and may be justified rigorously that in the 
second term on the right the order of integration may "be 
interchanged. Xhen with 



there follows as Integral equation for ^ , 



= »iT (106) 



lr(y*) « - JL.fjJL aTl* /TTl ^^(e.y*)a| (107) 

2 ^ 4sw^ dTl* y* - T|* J 1 - 1 '^o^i' .3^.'^^' 



-1 . -1 

Since, according to equation (86) 

= ir (108) 

it folloxirs that equation (107), the integral equation of 
lifting-line theory, determines directly the lift distrihu- 
tion in the stationary case. 
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In addition to the integral equation for the lift di8~ 
tribution further results may be deduced from equation (106), 
Introducing the value of uqi as given Tsy equation (105) 
into equation (87) for the section moment mQCx^i -1| l) 
there follows 

ao(xo. -1. 1) = »o(xo. -1. 1) - -^f Jirr^^'- ' **^'^=/pArfF 

(109) 

where mo stands for the value of the moment in the absence 
of an aerodynamic span effect. Vrom equation (109) follows 
that there is no aerodyi j iamio span effect for the moment 

about the quarter chord 'poinl^ » - ~ (which vanishes for 

2 

the straight line profile)- that is, 

ino(-|' -3-. 15 = mo(-|' 15 (no) 

A probably more important result of this nature concerns 
the value of the aileron hinge moment which is defined by 
equation (88). Introducing equation (105) into equation (88) 
gives 1 ^ 

mo(c.o.l) « mo(c.c.l) - Ax - c) /EZi dx f 

4Tre J vl + x Jy*«.T|* 

c -1 

Writing 

m n ' / .N J \ _i 

cos c 



f (c) -/ (x . c)y|ZI (i . |)AT7»- (1 * c) 



(112) 



and observing equations (107) and (108), there may be written 
instead of equation (ill) 

mo(o,c,l) - flo(c,c,l) ^ 1q " tp f (o) tp (lis) 
mo(c,©,l) Tp fio(e,o,l) 

Equation (113) indicates that t he aerodynamic span effect 
for the hinge moment differs from the span effect for the 
llf j> and in which manner the two are related. 



NACA TVS 2To. 946 



41 



She starting point is the complete integral eg^uation 
(89) v/hieh may be written in the form 

Wo « Ix + Is (114) 

where li stands for the integral over the airfoil region 

and la stands for the integral over the wake region. The 

appropriate approximation for I3. has heen obtained in the 

preceding section on the stationary airfoil. According to 
equation (lOl) 

I, . , i/ ^oCg .y) , l,/£.-i!L (U5) 

-"J ^ - i ^ 4TrJ d-n y - Tl 

Xt remains to obtain an approximation for Zg which 
corresponds to that obtained for I].. It will be shown that 
several essential steps are involved in the derivation of 
this approximation. 

Starting with the exact expression 

T _ ike^^ r r -ijctj r(n)(x - i ) 



X .s 



<|(x - 1)^ + (y - 



r'(Tl) / ^ > I . ,\ 

y - Ti[ ^(x - 1)« + (y - -n^'y^^ 



dldT] (116) 



the first step consists in separating from tg the value Ig'. 
would possess if the two-dimensional theory were correct. 
According to equation (94) the value of Is in the two-dimen> 
sional theory is 

- Ike^^-- - ".-ikl 
I3 



2tT 



/--IKS 
^ d| (117) 



Gi^he second step consists in writing 
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r(y) = If r.(Tl) 'y - dT) 



(118) 



which may readily Tae verified in view of the fact that 
r(±s) = 0. 

With equations (117) and (118) Ig of equation (116) 
may he written as 



Is = Is + 



r'(T|) 



ike 



eo s 

^ r -ike 



4rrt 



■J} 



1 MS 



- Tl /(x * I) + (y - Tl) 



r(ri)(x>- e) 

|(x . |)^ + (y - 71)^J^^ 

+ 1 - liLi-BL\t dUn (119) 



X - I 



i. further transformation is accomplished hy the follow- 
ing integration hy parte 

8 



r(Ti)d'n 



■n - 7 



|(x-|)^ i Z-r^fy^^ "^^'^ (x-|)V(x.|)^ 1 (y.^)^J"^ 
^ r'(?,)(7i 3 

J (x - trju - |)^ + 



y)d'n 
(y - r\) 



(120) 



The integrated part of equation (120) vani she b as T vanishes 
at hoth limits. Introducing equation (120) into equation 
(119), there follows 



ss Is + 



ike 



ik 



4tt 



1 *.e 



y - Tl 



y - Tl[j(x , t)^ + (y - n)^ \ 

X. J 



(x - e)/(x - + (y - Ti)*^J 

1 . iiLz4i \Ui 



— ^ y^ti6.y\ (121) 



and comhining the first two terms within the braces, 
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1-8 S ,S„- 



4Tr 

X - I 



The next step consists In separating in equation (122) 
the Integration with respect to C In two parts as follows: 



I 'I -I 



Q}hen 



where 



Ig = I3 + I4 + Ig (123) 



ike^^ r e-^^^ r(n) J , _ !y - tiI 



I4 = 



X -s 



+ - g) ^ (y - -n) ^ I ^ ^ ^ 

.-I J 



4Tr 



i! y - ^ - - I 



+ 



((x ~ |)^ 4- (y - 71)^} 



l/a-^ 



>d|cLT| (125) 

X - e ' 



Equation (124) may "be simplified by introducing a new vari- 
able of integration 



whence 

I4 



lc( I - x) = X, kdt = dX (126) 
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So far, no approximations have Tseen made in the treatment of 
the wake integral. The place where they are made is in the 
remaining term I5. Shis is possi'ble owing to the fact that 
the l-region of integration, (;c,l), in I5 is always small- 
er than the chord of the airfoil and that part of the factor 

multiplying e'^^'^T'Ct^) in the Integarahd is , as was the case 
previously in the integrals exitended over the airfoil region, 
of appreciahle magnitude in the nelghh.orhood of the line 

Tj = y only. Ihe terms in question are 

{•m/s 
(x - t)^ 4- (y . Tl)^) |y . Til ^ 1 X - j ^^^^^ 

(x - |)(y - 71) 2 (y - Tl)|y - n| 



and the *< sign holds as soon as (y Tll is somewhat larg- 
er than jx - || which, considering the assumption regarding 
aspect ratios as expressed by equatioi). (98) means over most 
of the span. With this oTiservatlon the approximation for I5 
is 



ike^^ p, f„\ r r e 



4tt 




In the first term of equation (129) the integration with re- 
spect to I may he carried out, the second term cancels he- 
cause the integrand Is an odd function of y - Tl. Hence 

.ik(l-x) ^ 



T ^ 1 - e^^^^-'^^ r r'(Tl) ,^ 

Ig « . 6 (130) 

4Tr J y - 



. y - 

— s 

Introducing I5 from equation (130), I4 from equa- 
tion (127) and I3 from equation (117) into equation (123) 
there follows 



ike^^ r e-^^t ^ ie^'^^^"^) / r<cn) 

Ztx J, X - I . 4tt J y - Tl 

r 1 "|i/8 -s 




1 + 




k|y - n| + k=(y - Tl)^f 1 1 - e^^^^-^'^ rr«(Tl) 
— L_ L _ , . ^ &.X>d.T\ + J ^ dTI 



^ J I 4Tr J y - T] 

-s 

(131) 
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Introducing Ig from equation (l3l) and Ii from eq^ua- 

txon (lis) into equation (114) there follows as approximate 
integral equation of lifting- surface theory for the rectangu- 
lar nonstationary airfoil 

1 00 

,ik ./^^-ikl 



-1 1 



- I 



4tt v4 y - n I ,^ L 

M /a 

- -1-^ ^ J^-M 



dT] (132) 



Putting as an abhreviation 



00 '. ( 



S(lc*y - Tl!) = i y^e-^^ [l ^ k I y - T| i , 1 ^ k^(y ^ Tl)^^ "j ^ ^ 



6^ 

equation (132) may be written in the form 



(133) 




00 



1 



2tt 

tim. A . S(kiy - .Tll)1 



S(kiy - Tl! ) j dT] (134) 



— s 



Equations (134) and (133) are the generalization of 
equation (94) of the t wo-rdimensi onal theory and of equation 
(lOl) of the three-dimensional stationary theory. They re- 
duce to these special cases when r'dl) = 0 or k = 0, 
respectively. It may "be noted that according to equation 
(134) the induced velocity due to the finite-span effect 
varies across the chord in contrast to the result of the 
stationary theory in which this velocity is uniform. In ad- 
dition to this the cumulative effect of the spanwise rate of 
change of r is modified as compared with the result of the 
stationary theory "by the occurrence of the function S which 
can "be tabulated once for all. 
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It is "because of the chordvlse variation, of the induced 
velocity that it is no longer -posei'ble to speak of a lifting 
line. To express the fact that the derivation of the inte- 
gral equation (lS4) depends on the assum-ption of the span "be - 
ing rather lon;;er than the chord the name of lifting- str it) 
theory is proposed for the basic equation (134) and for the 
consequences derived therefron . 

Further treatment of ecLuation (134) is possi"ble "by- 
means of a com"bination of the known procedures of the two- 
dimensional theory and of the three-dimensional stationary 
theory . 

As in the three-dimensional stationary theory an inte- 
gral equation for P is ohtained which may "be solved "by 
numerical methods. 

As in the two-dimensional theory explicit expressions 
for the* section lift and the section moments can he o"btained 
in terms of the function T. This will "be done in what fol- 
lows for lift and moment of the entire chord, leaving calcu- 
lation of the aileron hinge moment for future work. 

The first step in this program consists in the determi- 
nation of Uo from eq_uation (134) "by means of the inversion 
formulas (l02) and (103). The result is, if as a further 
a"bT5r eviation there is put 




(135) 



, eiiildl (136) 
4Tr J X - I 



The section lift and the section moments may "be calculated 
"by introducing equation (l36) into equations (86) to (88). 
"Wrom the two-dimensional theory it is known that these cal- 
culations lead to explicit results in terms of known func- 
tions and in terms of F for part of equation (136), It 
will "be shown that also the remaining terms can "be expressed 
in terms of known functions and in terms of (J. 
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2HB aqrjATIOH DBaJBRKIKISG THE OIROULATIQK JUNOTIOST 



Integrating equation (136) 
1 1 



lo evaluate eq.uation (13?) the order of integration with 
respect to x and | is interchanged and use is made of 
eq^uation (106). There follows 



It remains to evaluate the integrals 



CD 



-1 

1, = / / e"^^^ d| (140) 

-1 



■Im 



To oTjtaln Ig, interchange the order of integration 
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""^J?® °^ integral is known in terms of modified 

Bessel functions (See Durand, vol. II, p. 295) 

le = -Tr[^Ko(ik) + ICi(ik) - ±^ ] (143) 

?° o'btain I7 , use is made of the following known for- 
mula (see, for instance. Gray and Mathews Treatise on Bessel 
functions, p. 46) 

-ik| 



Jrom this there follows 



d| = TT Jo(k) (143) 



1 

= / -——1 e'^^^ dl =. ttF Jo(k) - iJi(k)l (144) 
-^1/1-1^ L J 

Substituting equations (142) and (144) in equation (138), 
1 

q, j^Jo(k) r iJi(k)'j (145) 



-1 

^ik 



and canceling the term on the left against one of the terms 
on tae right, rearranging and Introducing the value of 0 
from equation (135) there follows as integral equation for T, 

I riy) . . - . ^o(ic)-iJ,(k) / riM (1 . s(k, y - ,i ) ) d^ 

4ik[K:o(ik) + K;i(ik)3''i8 y - Tl ^ J 




.J V "F^.wod.y) H 

. - ^ - ^ (146) 

ike CKQ(ik) + K^(ik)3 
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A similar equation has 'been obtained "by Cicala (refer- 
ence 2*7) T3y means of considerations involving vortex fila- 
ments. Cicala's equation distinguishes itself from equation 
(159) "by the fact that the factor Jo(k) - iJi(k) is re- 
placed by a factor e~^^. Since for large k 

jo(k) - iJi(k) ^ /4- ^"^^ (^^'^^ 

V TTk 



it is seen that finite span corrections of a different order 
of magnitude in k are to be expected from Cicala's equa- 
tion and from the one given here. The difference is large 
even for relatively small values of k well within th© prac- 
tical range as can be seen by comparing the two terms as f ol- • 
lows, 

k = g sO.yS; e"^"*^ «0.7 - i 0;7; Jo(k) - IJi(k) = 0.85 " 1 O.36 (l^) 

It may be noted that Cicala's result would follow from the 
basic equation (134) if in the third term of this equation 
the factor exptik(l-x)] were missing. This means that in 
Cicala's work the chordwise variation of the finite-span- 
effect contribution to the induced velocity has been left 
out of consideration, lurther discrepancies which, it ap- 
pears, cannot all be accounted for in this manner are found 
between Cicala's expressions for section lift and moment and 
the expressions given in what follows. It may be emphasized 
that while there is a formal resemblance between part of the 
present results and Cicala's results, the present work and 
Cicala', s work are fundamentally different. Cicala's approach 
tct the problem does not permit a rational determination of 
all the factors of importance in the problem. 

Regarding KiSsBner'a work (reference 10) on the same sub- 
ject the following may be said. On the basis of the integral 
equation of lif t ing- sur f ace theory set up in terms of the ac- 
celeration potential Eussner obtains an integral equation for 
a quantity which« in the notation of the. present paper, is 
the section lift Ij^. Eussner states that this equation is 
correct only when Wo(x,y) = W{y)e~^^^ and for other func- 
tions Wq an equivalent w might be determined by means of 
the equation 
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It appears that Zussner' result cannot 'be correct inasmuch 
as it does not reduce to the value of the section lift of the 
two-dimensional theory when the assumption of two-dimension- 
ality is introduced into the equation. Turthermore, no pro- 
vision is made in Kussner's work for the determination of the 
finite-span effect on the values of the section moments. 

The comment that his result does not reduce to the ap- 
propriate two-dimensional result when it should do so applies 
also to Sears' special solution (reference 39) for the infi- 
nite periodically "bent airfoil. Sears' paper also contains 
no formula for the effect of three-dimensional flow on the 
values of the section moment. 

As Jones' work (references 31 and 32) deals with the 
transient problem of a rigid airfoil with elliptical form 
view no direct comparison is possible of his results with 
the ones given here. Inasmuch, however, as his developments 
make essential use of vortex-filament notions it appears de- 
siraTale to compare his results with results which may be 
obtained on the hasis of the notions of the present paper. 

DsicEHMiiirAirioisr of 7He bxfsession pgr ths sBC<rioiir lit^t 

According to eq.uation (86) 

Ife = iky^d - x) Uq dx + (149) 

and it remains to evaluate the integral 

1 

Ig = iky^d - x) uo dx (150) 
— 1 

Substituting eq.uation (136) in eq.uation (150) 

ik n -ik\ 




r r ik /-» -ik\ 

i_±J- IM— r r ±- — dX 

1 - I L ° 2-n J i - \ 



„ e^^(M> 1-^Ux (151) 



4'rT 
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Interchanging the order of integration with respect to I 
and X leads to the integral 



1 



-1 

Hence 



(152) 



-1 1 
^ikd^e ) 



4Tr 



^jdl (153) 



It is convenient (but not essential) to make use of equation 
(138) hy writing equation (163) in the form 

^ik(i-l) -| 

+ -2 d| (154) 

4tt J 

It remains to evaluate 

1 OO 

p ^ P p ^-ikX 

I9 = 



and 



yi y 1^ i?^] dl (155) 

1^0 =y yi - |3 e-iJ^^ d| (156) 



-1 

Ig is found hy interchanging the order of integration 



I9 



1 - -1 

00 



= -Tr/e -1 



dA 



NACA 9o. 946 52 



Integrating ty parts 

CO 



a/ 

Using eq[uatlon (l4l) and the known formula 



I, = . J 111! . ^ r .-i>=^ r 1 . /Eu , ) 



00 

^-ikX 



d\ = Eo(ik) (157) 
there follows 



The integral Iio is found "by means of the following 
fornula (see Bessel Functions Gray and Mathews, p. 46 - 
reference 35) 

J, (x) = . JL -1- / e"^""^ yi • |3 d| (159) 




.2 



-1 



where now T represents the classical T-funotion and 

r(|) '\J^ ■ 

Thus 1 

J, (k) 

Iio = / -y/i - 1^ e-^^£ de = TT — rr— (160) 



=y -/i ^ |3 e-ik^ dE = 



k 
-1 

Introducing now Ig from equation (158) and I^j from eq.ua- 

tion (160) into equation (154) there follows after some can- 
cellations 



Ig = - ik 



-/l - |a wod.y) d| 

-1 



iv-.ik „ r , , _-ik -) ivftils J,(k) 

^ — r Zi (ik) - -2 _ Q _3 (161) 

L ik J 4 k 



Substituting Ig from equation (l6l) in equation (149), 
there follows for the section lift 
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Ik = - ik J y 1 ~ |3 Wo(l.y) dl 

-1 

+ iM^K,(ik) r(y) - iiM (162) 

2 4 k 

where F is to "be determined from equation (146) and 0^ is 
defined Tsy equation. (135). Tor two-dimensional motion eq_ua- 
tion (162) reduces to 

1 

Ik = - Ik / yi - 1= VoCn 4l - Ki(ik) r (163) 

and, since according to equation (146) the two-dimensional 
circulation is given hy 




£ = - ^ ^ (164) 

2 1 Tr 

ike^'^CKo(ik) + Ki(ik)] 



there follows as special case the known expression for the 
section lift of the two-dimensional theory 

1 



J ° Ko(ik) + Ki(ik) ./y 1 - I 



-a -"^i 



(165) 



The factor of the second integral has been designated 
"by Theodorsen (reference 6) hy 0(k) = B'(k) + i &(k) . 

It is apparent that once the integral equation (146) 
for r has been solved the calculation of the lift distri- 
bution is no more complicated in the three-dimensional the- 
ory than it is in the two-dimeAsional theory. 



DBTEEMIirATIOEr OF THE EXPRESSION TOR THE SECTION MOMENT 



According to equation (87) the section moment about the 
semi-chord point is given by 
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1 1 

'^k^o^ ' j/^ (1 - x^) Uq dx +y^x Uq dx (167) 



-1 -1 

The second term of this eq^uatlon can "be expressed In terms of 
r and Ig of eq.uation (150) 

/ X uo dx. = - i2_ + L (168) 
-1 

Hence, vrlth equation (161). 

X Uo dx = / yi - |3 Wo dl - Ki(i]c)r + S- + ^ 

-1 -1 

4 k ■ 
It remains to evaluate ^ 

111 = ^ /"(l _ x«) Uo dx (170) 
-1 



/ 

-1 

which, according to equation (136) is given "b? 



•1 . -1 

00 



- J, / e ^ e ^ S. Idx (l7l) 

2tt J I X 4Tr J X - |J 

1 

Interchanging as "before the order of integration with regard 
to X and | the following integral has to he evaluated 
first 

1 



tTd - x^) /LZIi_i5_ = tt(|3 - I - (172) 
J 1 + X X - I ^ 3/ 



With this 
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00 




-1 



2/V 1 - e L Sit J I - X 



4tt 



] 



(173) 



WrUing 



and making use of eq^uation (138) In may be written in the 
form 



-1 



3tt 



4Tr 



(174) 



There remains to "be evaluated 

1 00 



-1 



(175) 
(176) 



The value of is found hy first interchanging the order 

of integration, 



dX (177) 



"by then integrating by parts 
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= IT - ^ 



OS 

-ik 



2ik 



(178) 



Utilizing equations (l42) and (15?) 

I 21k Ik.L Ik Ik ^ Ik ^1 ik J 

(179) 

Ihe value of Ij^^ is found with the help of equation 



(160) 

X 



-iBkJ dk ^ k ^ 

-1 



(180) 



Substituting now equations (179) and (180) in equation (174) 
there results after some cancellations 



+ Kidk) - ll ^ + Ko(ik)l + ^! Js(k) % (181) 

ik^ ik /J J 8 

Introducing equations (131) and (169) into equations (170) 
and (167), there followsj after some further cancellations 
and after making use of a recurrence formula for Bessel func- 
tions, as expression for the section moment distribution in 
the three-dimensional theory 

1 1 

= /TiT^ Wo d e - y^i yi - 1= Wo d i 

+ — — Ko(ik) T ^ - — J„(k) q. (182) 



mj5.(o) 



* ^ "0 

4 8 
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For the two-dimensional case this expression for the 
section moment reduces to 



1 1 



o _______ n _______ Ik ** 

£k(o) = Ul - |3 Wo dl - iiEj^/mr Wo dl + Ko(ik)r 

(183) 

and with F from ecLuation (164) 

1 r 



£k(o) =y yi - 1= Wo d| - ^J' \JT~^ w< 



-1 -i 

1 



. i Eoii^> / /XZI Wo d| (184) 

2 Ki(lk> + Kodk)^/^^/ 1-1 

This is in agreement with the known results for this Case. 

R^SUlSii or BASIC formulas op LirTING STEIP THEORY 
POR AIRJOILS OF REOTASTGULAR PLAIT lORM 

(l) The integral equation of lifting-surface theory as 
reduced for sufficiently large aspect ratio,, (s > 3) 

w,(..y) . _ 1 Ao^.^) -de . r^(y) ili^i dl 



s 



gikd-x) 



4tt 
where 



jT Ekli^l^l - S(kty - T|l)j dTl (134) 



1 



Vo = ikhk + h = e" bk, r = Tk e^*, i Tfe = / Uo dx 

dx >,3. J 

-1 

and 

en 



S(kly - nl) = i^y e-i^ [l ^ kly m , - A' ic^y - T) )^ 

0 



(133) 
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The function S is related to a function P which is tab- 
ulated In ecLuation (27), in the following way 

S(x) = ix P(x) 

(s) The integral representation for the chordwise veloc- 
ity component at the airfoil 

.ik n -ikX 

uoU.y; = ^ / - - d / " • > <wa - ^ — Tir J - — - a^ 



8 

—8 

(3) The integral equation for the circulation function 
3 ^ 4ikCKo(ik) + Ki(ik)3 ./ y - T) L J 



^/pX wodi 

^ ■ ^ (146) 

ike^^CKo(ik) + Ki(ik)] 



The modified Beesel functions of the second kind may "be ex- 
pressed in terms of Hankel functions,* 

Ko(ik) = - i 5 Ho^*^(k). Ki(ik) = - Hi( = )(k) 



(4) The formula for the section lift I = Ij. e 



it 



Ij^ = - ik / V 1 - &3 d| + K^(ik) Ty^iy) 

'^-^ s . ' 

. ikei^ / r^lM (1 , si dTl (162) 

4 k J y , r, \ J 
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(5) The formula for the section moment alsout the semi- 
chord point m(o) = mjj(o) e^^ 

mjj.(o) = y V^l - |3 dl - ^ - |3 dl 

4. Eodk) ric(y) 1- li^ Jo(k) <r ^^^l^(l - s] dTl (182) 

4 8 J y - Tl L J 



-s 2^ 



ooi)rcLUsiiir& bs marks 



The developments of this paper indicate that use of the 
velocity potential is preferalile in important respects to 
the use of the acceleration potential in thin-airfoil theory. 

Becognitlon of this fact is considered as an essential 
aid in the e st ahli shment of "lifting-strip" theory for a non- 
stationary airfoil. While this theory has here been devel- 
oped for airfoils of rectangular plan form it can 1oe extended 
to tapered airfoils. 

Application of the results of the present paper should 
permit, among other things, the investigation of the aero- 
dynamic span effect in the problem of wing flutter. Tor this 
purpose it remains to establish a convenient scheme for the 
numerical solution of the integral eq^uation for the circula- 
tion function. 



While some investigators have stated as their opinion 
that the error in flutter calculations resulting from the 
assumption of two-dimensional flow is neglibibly small for 
all wings with aspect ratios above three, the author consid- 
ers the available evidence as inconclusive. He believes that 
a decision on this oLuestlon can be reached by applying the 
results of part III of this paper to a number of representa- 
tive flutter cases. 

It is emphasized that the manner in which lifting-strip 

theory is obtained indicates that its range of validity is 
no less than the range of validity of lifting-line theory 
for the stationary airfoil. Inasmuch as experiments have 
shown that lifting-line theory may be applied for wings with 
aspect ratios as low as three, the same must be true for the 
results obtained in part III of this paper for, the nonsta- 
tionary. airfoil . It should be possible to apply lifting- 
strip theory to tail flutter problems. 

Massachusetts Institute of Technology, 
Cambridge, Mass., Teb. 1943. 



NAOA Til ITo, 946 



60 



1, Prandtl, L, ; Tragf lugeltheor ie» I, Mittellung, 

Kachrichten der K, G-esellschaf t der Wissenschaf ten 
zu (Jo'ttingen, MatJi,— phys , Klasse, 1918, pp. 451—477, 

2, Munk, Max M. : G-eneral Theory of Thin Wing Sections, 

Rep. No. 148, NACA, 1922, 

3, Birnbaum, W, : Das ebene Problem des schlagenden ?lugels. 

Z.f,a,M,M. , Bd. 4, Heft 4, Aug, 1924, pp. 277-292. 

4, Wagner, Her'bertl liber die Sntstehung des dynamischen 

Auftriebes voii Tragf lugeln. Z.f.a.M.M. , Bd. 5, 
Heft 1, Pe"b. 1925, pp. 17-35. 

5, &lauert, H, ; . The force and Moment on an Oscillating 

Aerofoil. R, & M. No. 1242, British A.fi.O., 1929, 

6, Eheodorsen, Theodores General Theory of Aerodynamic 

Instability and the Mechanism of Flutter. Rep. No. 
496, NACA, 1935. 

7, Cicala, Placido: lie Azioni aerodinamiche sui profili 

di ala oscillanti in presenza di corrente uniforme. 
Mem, R, Accad* Sci, Torino, ser, 2, pt, I, t, 68, 
1934-35, pp. 73-98. 

8, Ellenberger , G-, ; Bestimmung der Luftkrafte auf einen 

ebenen Tragflugel mit Querruder. Z.f.a.M.M,, Bd, 16, 
Heft 4, Aug. 1936,' pp. 199-226, 

9, Borbely, v.: Mathemat ischer Beitrag zur Theorie der 

Plilgelschwingungen, Z.f.a.M.M,, Bd. 16, Heft 1, 
ffeb, 1936, pp. 1-4, 

10, Edssner, H. ; Zusammenf as sender Bericht uber den 

instationaren Auftrieb von Plugeln, Luftfahrt— 
forschung, Bd, 13, Nr. 12, Bee, 30, 1936, pp. 410-424, 

11, G-arrick, I, B. s Propulsion of a Flapping and Oscillating 

Airfoil. Rep, No, 567, NACA, 1936. 

12, von Karman, Th, , and Sears, W, R, ; Airfoil Theory for 

Non-^Unif orm Motion, Jour, Aero, Sci,, vol, 5, no, 10, 
Aug, 1933, pp, 379^390. 



HAG A ITo, 946 



61 



13, Dietze, : Die Luftkrafte der harmonisoh schwingeudeix, 

in aich verf ormtaren Platte (Etenes Problem), Luft— 
fahrtf or ©Chung, Bd, 16, Lfg, 2, Peb» 20, 1939, pp, 
84-96, 

14, Schwarz, I|, i Berechnung der Druckverteilung einer 

harmonisch sioh. verf ormenden Tragflache in. ebenor 
Stromuag. Luf tf ahrtf orschuug, Bd, 17, Iifg, 11/12, 
Dec. 10, 1940, pp. 379-386. 

15, StShngen, Heinz; Bestimmung der Auf triebaverteilung 

fur beliebige Instationare Bevregungen (Sbenes 
Problem). Luf tf ahrtf orschung, Bd. 17, Ifg, ll/l2, 
Deo, 10, 1940, pp, 401-420. 

16, Possio, 0,: Ii'Azionc Aarodluajnica sul Profile 

Oecilla&te in un Pluido Compressible a 7eloclta 
Iposonora, I'Aerotechnica, vol, XVIII, faso, 4, 
1938, pp, 441-458, . . 

17, Blenk, H. » Der Blndecfcer als tragende Wlrbelf l."che, 

Z,f,a,M.K,, Bd, 5, Heft 1, Peb, 1925, pp. 36-47, 

18, von Karzuan, Th,, and Burgers, J, I-i^ t C-eneral Aero- 

dynamic Theory — Perfect Pluids, Vol,. II, div, E, 
eh, II of Aerodynamic Theory, P, Durand, ed, , 
Julius Springer (Berlin), 1935, pp„ 3-4-99, 



19, von Karman, Ih, j Keue Darstellun/; der Tragf lugeltheor ie, 

Z.f.a.K.K, , Bd, 15, Heft l/2, Peb, 1935, pp. 56-61, 

20, Schlichting, H, 8 Airfoil Theory at Supersonic Speed, 

T«K. Ko, 897, KAOA, 1939, 

21, PUchs, R, J Heue Behandlung der Tragf lugeltheor ie, 

Ing.-Archiv, , Bd. X, Heft 1, Peb, 1939, pp. 48-63, 

22, Bollay, tfilliam: A Non-Linear Wing Theory and its 

Application to Eectangular Vfings of Small Asyoct 
Ratio. Z.f.a.K.K, , Bd. 19, Nr. 1, Peb. 1939^ pp. 
21— 36 , 

23, Wieghardt, Karl: Ohordwise Load Distribution of a Siiaple 

Rectangular Wing, T,M. 1*0, 963, lUOA, 1940. 

24, Sinner, vr. J Die kreisformige Tragflache auf potont ial- 

theoret ischer Grundlage, Ing.-Archiv, Bd, 7II1, 
Heft 1, Peb, 1937, pp, 47-80, 



KAOA EST Ko, 946 



25, Krienes, Klausi The Elllptio Ving Based on th. 

Theory. T.M. No, 971, iTACA, 1941. 

26, Cicala, P.: Sul moto uon stazionario di un'ala di" 

allungamento finite, Bendicontl R, Accad, Naz, 
Lincei, vol, 25, ser, 6a, no. 3~4, ' Aug, 1937, 
pp, 97-102. 

27, Cicala, P.: Comparison of Theory with Experiment " in 

the Phenomenon of Wing Flutter, T.M, ETo, 887, 
IIAOA, .1939, 

28, Borbely, v. j Uber einen Grenzfall der Instat ionaren 

raumlichen Iragf lugelstromung, Z.f.a.M.M^, Sd, 18, 
Heft 6, Bee, 1938', pp. 319-342, 

29, Sears, William s A OontriTjut ion to the Airfoil Theory 

for Sonr-TJnif orm Motion, Proc, Fifth Int. Cong, Appl, 
Kech,, 1938, pp, 48Jt*487, 

30, Kussnor< H. a.: General Airfoil Theory, T,K, Ho, 979, 

IIACA, 1941, 

31, Jones, Robert T, : The Unsteady lift of a Tinita Wing, 

T.N, No, 682, NAOA, 1939, 

32» Jones, Roljert T, ? The TJnsteady Lift of a Wing of Finite 
Aspect Ratio. Rep, No. 681, NACA, 1940. 

33, Prandtl, L< S Recent Work on Airfoil Theory. T,M, llo, 

962, NACA, 1940. 

34, Sohngen, H. } Bie Losungen der Zntegralgleichung 

g(x) ss -i. r — ~-. d| und deren Anwendung in der 

Tragf lugeltheorie. Math. Z, , vol, 45, 1939, pp. 
245-264, 

35, aray, A., Mathews, 9, , and Macfohert, T. M. : A 

Treatise on Bessel Functions and Their Applications 
to Physics. MacKillan and Co., Ltd., 2d ed, , ISSl^ 



LTAOA !?iT STo^ 946 



63 



BIBZ.Z Oa&AFKY. 



Eildebrsnd , 3« , and Orout, P, D, : A Least Square 

Procedtt-re for Solring Integral Equati-ons by Polynomir.l 
Approximations. Jour. Math, Piiys., vol, 20^ 1941, 

Hildebrand Prancis B.J A Least-Squares Procedure for 
the Solution of the Lif t ing~Iiine Integral Eq_uation, 
TIT "0, 925, 1944, 

Blot, i:. A,: Some Simplified Methods in Airfoil Theory, 

Jour, Aero, Sci,, vol. 9, no, 5, March 1942, pp, 185-19C, 

Oicala, P,! Zuschriften an den Herausgeber, Z,f,a,M,i:, 
Bd, 19, Nr. 6, Dec» 1939, p, 3S4, 

Oioala, P,J Zuschriften an den Eerausgeber, Z, f , a, M, i-i, , 
BA, 20, iTr. 3, June 1940, p, 184, 

Xussner, H, t Zuschriften an den Herausgeber, Z, f , a, K, I-.', , 
rsd. 19, Nr, 6, Deo. 1959, p, 334^ 

Kvissner, H, &, s Zuschriften an den Herausgeber, Z, f ,a« IT, , 
3d. 20, iTr, 3, June 1940, p. 184. 

Schade, T'h. J Theory of the Oscillating Circular Arc Air- 
foil on the Basis of the Potential Theory, Part I, 
(Translation^), Luf tf ahrtf orschung, 3d. 17, Lfg. ll/lS, 
Doc, 10, 1940, pp, 387-400, 



Available for reference or lonn in the Office of 
Aer onau-o leal Intelligence, ITAOa, Washington, C, 



"SXCA, Technical Note Ho. 946 




